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INTRODUCTION 


There is a wealth of experimental data on the transport properties of gas 
mixtures - diffusion coefficients, mixture viscosities and heat conductivities, 
as well as results on thermal diffusion. Rather then attempting some sort of 
statistical assessment of all this information, a few simple systems which 
typify interactions amongst various types of molecules - monatomic and poly- 
atomic, nonpolar and polar - are examined in the hope of obtaining some insight 
as to the concordance between theory (rigorous or approximate) and experiment 
(painstaking or prefunctory) . 

To this end, five binary gas systems are considered in the neighborhood of room 
temperature, since precise data are most easily obtained here. For the systems 
chosen, there is data on the three most important transport properties: the 

binary diffusion coefficient, mixture viscosity, and thermal conductivity. 
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There is also thermal diffusion data for some of these mixtures, but it is not 
analyzed or considered in detail here. 

Systems have been selected so as to avoid undue complexity - ternary or multi- 
component mixtures are not examined, since the behavior of these systems should 
be understandable in terms of the behavior of binary systems. Similarly, 
systems involving interacting pairs of complex polyatomic (and possibly polar) 
molecules have been avoided. 


Sn general we expect that the rigorous Chapman-Enskog theory for monatomic 
gases should serve rather well in describing the diffusivity and viscous prop- 
erties of polyatomic gases. In the case of the thermal conductivity, however, 
a diffusive flux of internal energy (rotation, vibration, etc.) can make a sub- 
stantial contribution to the heat flux, and must be considered. To this end 
the recent approximate theory of Mason and Monchick' for the heat conductivity 
of polyatomic and polar gases is invoked and extended (in an approximate 
fashion) to gas mixtures. 

Specifically, the following mixtures have been examined: 



Hel ium-Arqon . Since this is a mixture of monatomic gases, obeying a 
spherically symmetric intermolecular force law, we should expect the 
rigorous Chapman-Enskog treatment to be fully applicable. 

Hel i um-Hydroqen . A mixture of a monatomic gas with a diatomic molecule 
whose only internal energy is rotational. The transfer of energy be- 
tween translation and rotation is slow for hydrogen, and requires 
several hundred collisions. Thus we should be able to treat the trans- 
port of translational and internal energies separately. 

Hydroqen-Ni troqen . A mixture of diatomic gases which possess rotational 
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internal energy only. However, in examining the heat conductivity, 
the translation-rotation relaxation in nitrogen must be considered. 
Argon-Ammon ia . This mixture involves the interaction of a monatomic 
gas with a gas possessing an appreciable dipole moment. 

Hyd rogen-Ammon ? a . A mixture involving a polar molecule and a diatomic 
molecule. 

In general, the emphasis here is on mixture properties; hence we try to tie to 
perfect predictions as to the properties of the pure components. To this end 
experimental thermal conductivities were corrected to be in accord with theory 
for the pure gases for those substances where the theory seems sound and near 
rigorous - helium, argon, and hydrogen. Experimental conductivities of nitro- 
gen and ammonia have been assumed correct, and a rotational relaxation time 
has been assigned to make the experimental heat conductivity of nitrogen or am- 
monia correspond to the theory of Mason and Monchick*. 

THEORETICAL CONSIDERATIONS 

In general, the transport properties of dilute gas mixtures involve terms char- 
acteristic of the interaction of the pure components with themselves, and also 
all possible pairwise interactions between unlike species. The self-interactions 
give rise to the properties of the pure species - not the subject of this dis- 
course. The parameterwh i ch most directly characterizes the unlike interactions 
is the 

Diffusion Coefficient . Such data have been analyzed in terms of the first 

Chapman-Enskog approximation (Reference 2, Equations 8.2 - 44), in most 
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cases corrected by an average of the Kihara second approximations cal- 
culated for the extremes of composition. 


4 


Mixture Viscosities have been analyzed in terms of the first Chapman-Enskog 
approximation (Reference 2, Equations 8.2 - 22). However, in most cases the 
input quantity ^ ^ derived from diffusion has been corrected according to 
the third C-E approximation for pure gases . 

Mixture Thermal Conductivities for monatomic gas mixtures have been obtained 
from the expression of Muckenfuss and Curtiss**. This expression was also used 
to compute the translational energy contribution to the conductivity for the 
other mixtures, but with an appropriate correction, as discussed below. The 
internal energy for mixtures has been computed according to H i rschfelder 1 s'* 
Eucken-type approximation for mixtures. Again with an appropriate correction. 
The use of the Muckenfuss-Curtiss equation for the translational heat conduc- 
tivity with Hirschfelder's equation for the internal heat conductivity is 
justified assuming inelastic collisions rare. This assumption is necessary 
in order that the translational distribution function should not be unduly 
perturbed so that the translational conductivity is related to the viscosity 
or diffusion coefficients as in the case of the noble gases. It seems that 
this assumption is justified, except in the case of rotation-translation in- 
terchange which occurs every few collisions (say 2 - 20). 

Mason and Monchick' have derived approximate expressions for the heat con- 
ductivity of pure polyatomic gases taking account of inelastic collisions. 

They find that the translational conductivity is less than that of a monatomic 
gas: 
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Here C trans , C rot are the translational and rotational portions of the heat 
capacity while Z is the collision number for rotational relaxation. 

On the other hand, the internal energy contribution to the conductivity is 
enhanced: 



, a /s 1 

Tr\ x ijCrf Z 


(2) 


As yet no comparable theory for mixtures has been published. Consequently 
we use an extension to binary mixtures arrived at intuitively as follows: 
Assume that the mixture conductivity may be written 




(3) 


where , A are the deviations attributable to the separate components, 

An approximate formula^ for the conductivity of a monatomic gas mixture is 


'Xw, 






I > Q’Voyw 


. "u r « 


r yp-) r M9 , K 
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where are the conductivities of the pure components, , 9^ 

are mole fractions, and are functions of the molecular weights and 

cross sections^, H i rschfelder 1 s formula for the internal heat conductivity is 
of the same form as Equation (4), but with ^,4", Modified Eucken in place of the 


^•mon_ an d the N'T, replaced by Di/Di j. We assume that in a mixture the terms 2' 
0 t 


-1 


must be replaced by 
-1 




-1 


( 5 ) 
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Hence, from Equations (1), (2), (4), (5) and Hi rschfelder's internal conductivity 
formula, we infer 


A/ll = ~(£ rot 


h + k 
2, 


^ I j 1*^ MoUrfi ec L / Cl,. 4 A i ^OYloK. /0<j. 


1 


(6) 


_l__ D i 1C- a_ 
L- 1 Dia-Xl 


- ^ 


The formula for A^is obtained from Equation ( 6 ) by interchanging subscripts 
1 and 2. Thus Equations (3) and ( 6 ) were used in analyzing data on mixtures 
involving nitrogen and ammonia. 

As mentioned in the introduction, mixture conductivities were corrected so as 
to match theoretical values for helium, argon, and hydrogen. Hot wire thermal 
conductivity cells measure in essence the reciprocal of the conductivity. Thus 
if A | , A ^ are the experimentally reported conductivities of the pure com- 
ponents, and /^'m,is the experimental mixture value, while ^ ( , "A^are the 
theoretical conductivities of the pure components, then the corrected mixture 
conductivity is 


n “ * __ / ^ "‘‘l \ ^ 3 . /I I 

/ou - L A l ) ^ 


-1 


+- a 


-1 


(7) 


RESULTS AND DISCUSSION 

Hel i um-Arqon . Diffusion coefficients for the helium-argon system were calculated 
assuming the exponential -6 potential with the force constants = 3 . 505 ^. 

06 - 13.21, and 6 A. - 33.4°K. The values for 06 and e/k. are taken from Mason^ 
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while V^was adjusted to fit the experimental data of Walker and Westenberg^ 
at 25°C and 1 atmosphere. Comparison between theory and experiment are shown 
below: 

TABLE ± 

( 2 ) 

Helium Concentration Diffusion Coefficient f D 

Experiment Calculated 

0.0 (trace He) 0.754 0.750 1.0284 

1.0 (Trace Argon) 0.725 0.729 1.0001 

( 2 ) 

The factor f D is the ratio of the second Kihara approximation to the first 
Chapman-Enskog approximation. Of course the agreement between experiment and 
theory is not surprising, since Y^,was chosen to obtain a good fit. The ex- 
perimentally observed variation of diffusion coefficient with concentration is 
somewhat larger than that predicted theoretically, but it is not clear if this 
is a real effect or merely experimental error. 

Viscosities for helium-argon mixtures were calculated again assuming the 
(exp-6) potential with force constants as follows: 

TABLE 2 

w 06 e/k. 

3.109 12.4 9.16 

3.858 14.0 123.2 

3.471 13.21 33.4 

3.494 " " 


He - He 
Ar - Ar 

He - Ar (combining rules) 
(Diffus ion) 
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Again the and €^lc are due to Mason' 5 ; for pure helium and argon lf^_has been 
chosen to fit the recent, precise viscosity determinations of Iwasaki and 

p 

Kestin . Two values of ^ for the unlike interaction have been tested - the 
first from Mason's^ empirical combining rules, and the second obtained from 
the experimental diffusion coefficients as follows: 


^ ViScoSiti 


a 


% 

= 





( 8 ) 


(3) 

Here f ^ is the correction for the third approximation for the viscosity of 
a pure gas (in this case 1.0062). This procedure lacks elegance, but is prob- 
ably in the right direction. Theory and experiment are compared in Figure 1. 
Although the combining rule is not bad - the worst deviation is only 0.7% - the 
force constants derived from diffusion are clearly to be preferred. All in all, 
the agreement between the diffusion and viscosity data is heartening indeed, 
and seems to fully confirm the rigorous Chapman-Enskog theory. In fact the 
viscosity data are apparently of a quality such as to justify a proper calculation 
of the second approximation. 

Thermal conductivities have been computed using the force constants of TABLE 2, 

( = 3. ^9^ for the He-Ar interaction); deviations of experiment from theory 

are shown in Figure 2. Clearly the deviations are much more serious, amounting 
to 1.5% in the best case (the data of Von Ubisch), with errors approaching 10% 
at worst (Baker and Thornton). In an overall sense, the errors seem random 
rather than systematic, and there seems no cause for questioning the rigorous 
theory. Indeed, the computed conductivities are perhaps more reliable than the 
experimental values for this system. 
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ii 

Data on thermal diffusion in helium-argon mixtures has been analyzed by Mason ; 
agreement between theory and experiment seems satisfactory. 

Hei ium-Hydrogen . Diffusion coefficients have been measured by Bunde^® with 
what appears to be good precision and accuracy. Here they have been calculated 
assuming the exponential -6 potential with force constants .289^, 

OC = 13. 22 , and €/|c = 18.27°K ( and s/i. again from Mason^), with (^adjusted 
to bracket Bunde's data at 25° and 1 atmosphere, as shown below: 



TABLE 2 



Diffusion Coefficient 

f^D 

Range of 

Helium Concentration 

Experiment Calculated 


0 (trace He) 

1.543 

1.0036 

0 - 1.00 

1.549 

- - - 

0.1490 - 1.000 

1.550 

- - - 

0.4537 - 1.000 

1.572 

- - - 

1 . 00 (trace H 2 ) 

1.571 

1.0221 

Again the experimental 

variation of diffusion coefficient with 

compos i t i on seems 

somewhat larger than theory, but in this case the experimental 

compositions em- 

brace a range of values 

, so that the comparison is less meaningful. 

Viscosities and thermal 

conductivities for helium-hydrogen were calculated as- 


suming force constants as follows: 
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TABLE 4 


I ow 


06 


A 


He - 

He 

3.100 

12.4 

9.16 

H 2 " 

h 2 

3.337 

14.0 

37.3 

He - 

H 2 (combining rules) 

3.226 

13.22 

18.27 


(0 i f f us i on) 

3.278 

1 1 

II 

Again oi 

and *4- are from Mason, 

with for the 

pure components 

adjuste 


as to reasonably reproduce the experimental data for these gases, and for 
the He - H 2 interaction has been adjusted according to Equation (8) with 
f^^ = 1.0063. Theoretical viscosity calculations are compared with experiment" 

in Figure 3. Again the agreement between the diffusion and viscosity data is en- 
couraging. (We should note, however, that the viscosities seem considerably less 
precise than in the helium-argon case, with deviations of 1 - 2% rather than 0.2%) 
Thermal conductivities computed according to the modified Eucken approximation for 

I A 

mixtures are compared with Barua's z experimental results in Figure 4. (The modi- 
fied Eucken approximation assumes that the collision numbers and Z^-He are 
sufficiently large that the correction of Equation (6) is negligible.) Clearly 
these data seem to favor the empirical combining rules over unlike interaction 
force constants determined from diffusion coefficients or mixture viscosities. 

(it might be noted that these experimental data are in even closer accord with 
the theoretical analysis accompanying the data.) 

However, from research on gas chromatography using thermal conductivity detectors 
it is well-known that small amounts of hydrogen decrease the heat conductivity of 
helium-hydrogen mixtures . From this fact alone one can establish a maximum 
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slope to the conductivity-concentration plotatzero hydrogen concentration; 
these slopes are shown in Figure 4. Further, the data of Schmauch and 
Dinerstein^ indicate that a mixture containing 17.3% hydrogen has the same 
thermal conductivity as pure helium; the consequences of this fact are in- 
dicated as an "X" on Figure 4. 

Summarizing, it seems that the rigorous Chapman-Enskog theory properly repre- 
sents the relation between diffusion coefficients and viscosities of helium- 
argon mixtures. It is not clear whether or not the modified Eucken approxi- 
mation describes the thermal conductivities of these mixtures. Experimental 
data of enhanced accuracy are clearly required. 

Thermal diffusion data have been analyzed in reference 2., p 584. 

Hydroqen-Ni trogen . Diffusion coefficients have been measured by Bunde'® using 
a Loschmidt-type cell; he obtained no information on concentration dependence. 
Values have been calculated assuming W = 3.686^ ( to fit the experiment), 

06 - 15.56 ande/Jc = 58.97°K (f rom Mason^). Comparisons at 1 atmosphere are 
as follows: 

TABLE £ 

Hydrogen Concentration 
0 0-1 1 


Diffusion Coefficient 

(calculate) 

(Experiment) 

(calculated) 

25°C 

0.7842 

0.7835 

0.7599 

55° 

0.9231 

0.9079 

0.8944 

85° 

1.071 

1.052 

1.037 


f (2) D (25°) 


1.03220 


1.00006 
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The correction for the second approximation was calculated at 25°C and then 
applied at the other temperatures as well - however, this factor should vary 
only very slightly with temperature. 

Viscosities and thermal conductivities for hydrogen-nitrogen mixtures were 
calculated from the following force constants: 


TABLE 6 

06 */k_ 


h 2 - h 2 

3.337 

14.0 

37.3 

n 2 - n 2 

4.011 

17.0 

101.2 

H 2 - N 2 (combining rules) 

3.690 

15.56 

58.97 

H 2 - N 2 (Diffusion) 

3.675 

1 1 

II 


All values were taken from Mason^ except for V^from diffusion, again adjusted 
according to Equation (8) with f^ = I.OO 63 . In Figure 5 experimental vis- 
cosities^*^ are compared with theoretical values calculated us i ng 
from diffusion. (Values from combining rules were only slightly different 
and hence are not shown.) The agreement between viscosity and diffusion data 
again seems satisfactory. 

Thermal conductivities have been computed from Equations (3) and (6) with 

Z u = z ,. 1 , = 00 and Z., = Z„ u =7.2. Theory and experiment ^ are com- 

H 2 H2”N2 N 2 N2“ m 2 

pared in Figure 5. Also shown as a dashed line is the modified Eucken-type 
approximation corresponding to Z^ = Z^ ^ ~ 00 * * t wou ^ seem that the 

apparent discord between theory cannot be rationalized by the choice of larger 
collision numbers. Once again better experimental conductivity data would be 


most welcome. 
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Argon-Ammonia . Diffusion coefficients for this system have recently been de- 
termined by Srivastava and Srivastava^. They fitted their experimental data 
to the Lennard-Jones (12-6) potential withd” = 3.286$ and^/k = 224.65. These 
values agree quite closely with estimates based on empirical combining rules 
using viscosity force constants of the pure components ( <T = 3.271$ and^/k 81 221°K). 

Viscosities of argon-ammonia mixtures have been measured by Iwasaki, Kestin, and 
1 8 

Nagashima . These data are analyzed here by choosing viscosities of pure am- 
monia and argon and calculating the parameters characterizing the unlike in- 
teraction from the force constants derived from diffusion; results are shown in 
Figure 6. Clearly there is a serious discrepancy, with experimental viscosities 

as much as 7% larger than theory. It is possible to obtain agreement by arbi- 

o 

trarily takingO"* = 3.072A for the unlike interaction. This is equivalent to 

assuming that the measured diffusion coefficients are too low (should be some 

14.4% larger). The lower portion of Figure 6 shows that accord between theory 

and experiment is much improved; errors are generally in the range of 1 - 2%. 

(it might be noted that these experimental data appear much less precise than 

8 

the same authors' results for the helium-argon system.) 

In view of the serious discrepancy between the viscosity and diffusivity data 
we set about to measure the heat conductivities of these mixtures at 300°K. 

A hot-wire thermal conductivity was used and calibrated with helium and argon 
assuming conductivities of these gases to be 370.9 and 42.42 microcalories 
cm” ' sec” ^ °K” ' respectively. (These values were obtained by considering liter- 
ature data on the thermal conductivity and, more especially, viscosity of pure 
helium and argon .) Data were analyzed using Equations (3) and (6) after 
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correcting the coefficient for the diffusion of internal energy for the resonant 

exchange of rotational energy postulated by Mason and Monchick*. The collision 

number Z. 1U was taken as 2.26 so as to fit the pure ammonia datum. An upper 
NH j 

limit to the mixture conductivities was calculated assuming - Ar = CO . 

Calculations are compared with experiment in Figure 7 - it is seen that the ex- 
perimental conductivities are compatible with the viscosity results but incom- 
patible with the diffusion coefficients. As a matter of fact by assuming 

Z NH 3 “Ar = * ,ifZ NH 3 the ex P erimental data can fitted almost perfectly. 

Thus the experimental diffusion coefficients*^ for the ammonia-argon system 
are clearly suspect and should be remeasured. 

Hydrogen-Ammonia . Diffusion coefficients for this system have been measured 
by Bunde ; his results have been fit assuming a Lennard-Jones (12-6) potential 
with C3~= 2.933A and = 142.7. (The value of e/fc. was obtained by Bunde from 
a fit of the first Chapman-Enskog approximation to his data over a range of tem- 
perature; <3~ has been deduced by applying an average correction for the second 
Kihara approximation to the diffusion coefficient to (T' - 2.927 deduced by 
Bunde.) Comparisons at 1 atmosphere are as follows; 


Diffusion Coefficient 


TABLE 2 


Hydrogen Concentration 
_ 0 _ 0 - 1 1 
(calculate) (Experiment) (calculated) 


25°C 

0.7891 

0.7830 

0.7825 

55° 

0.9403 

0.9426 

0.9324 

o 

LA 

OO 

1.1018 

1.0933 

1 . 0926 


f (2) 0 (25°) 


1 . 00846 


1.00005 
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The correction factor calculated at 25°C was applied at all three tem- 

peratures. 

Viscosities and thermal conductivities for hydrogen-ammonia have been analyzed 
by assuming viscosity cross sections to exactly fit the viscosities of pure 
hydrogen and pure ammonia. Parameters characterizing the unlike interaction 
were taken from the combining rules of reference 2, p 600 ( CT = 3.022, 

G /k.= 123°K) and also from diffusion ( <T - 2.931,*/^ = 142.7). (The value of 

0~~ from diffusion has been corrected according to Equation 8 with fl^ = 1.0015.) 

20 

Experimental viscosities of Trautz and Heberling are compared with theory in 
Figure 8. The agreement is seen to be very good, with the diffusion force con- 
stants somewhat to be preferred. Thus it appears that the rigorous Chapman- 
Euskog theory for monatomic gases is quite acceptable even for mixtures involving 
polar gases. 

Thermal conductivities were calculated from Equations (3) and (6) with the resonant 
correction of Mason and Monchick and assuming Z|^ =2.11 so as to fit the experi- 

mental value of Gray and Wright'^ for pure ammonia at 25.3°C. A theoretical up- 
per limit to the heat conductivities was calculated assuming = OO . Cal- 

culations are compared with experiment in Figure 8. Also shown as a dashed line 
is the modified Eucken-type approximation corresponding to =00. In this 

case the correction for inelastic collisions it quite large, and amounts to about 
15% for pure ammonia. It does not seem possible to draw any firm conclusions 
concerning the theory from Figure 8. The deviations are certainly larger than 
one would hope, but show a large random, rather than systematic, fluctuation. 

The agreement is scarcely worse that in the case of the H£ - N£ data of the same 
authors shown in Figure 5. 
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SUMMARY AND CONCLUSIONS 

Experimental transport properties of five selected binary systems have been 
investigated in terms of the best theoretical formulations presently avail- 
able. This analysis leads to the following tentative conclusions: 

Regarding Theory : 

1. The diffusion coefficients and viscosities of gas mixtures, 
including mixtures with polyatomic and polar molecules can 
be very well described in terms of the Chapman-Enskog theory, 
despite the fact that this theory is strictly rigorous only 
for monatomic gases. 

2. There is no reason to doubt the validity of the Chapman-Enskog 
theory for the thermal conductivity of mixtures of monatomic 
gases; however, none of the experimental data are good 
enough to provide a really definitive test of the theory. 

3. There are approximate methods which take account of the effects 
of internal energy on the heat conductivity of mixtures involv- 
ing polyatomic and polar gases. These formulations seem prom- 
ising, but again, experimental precision is too low to permit 

a meaningful evaluation of these methods. 

Regarding Experimental Data : 

1. There is already a large amount of mixture data of modest ac- 
curacy and precision (say errors larger than 1 - 2%). Further 
experimentation of this sort will do little to further our testing 
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and understanding of theory. 

2. The very best determinations of diffusion coefficients and mixture 
viscosities seem good enough to provide really meaningful tests 

of theory, including the second Chapman-Euskog approximation for 
diffusion and perhaps viscosity as well. More data of this 
quality would be most welcome. 

3. There is a crying need for thermal conductivity data of high ac- 
curacy and precision - little or none exists. It is proposed 
that such data may be acquired from precise measurements relative 
to the noble gases. Accurate conductivities for the monatomic 
gases can be computed from the best viscosity measurements by 
means of rigorous kinetic theory. 
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